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Abstract 

If the Boltzmann-Gibbs state ujn of a mean-field iV-particle system 
with Hamiltonian verifies the condition 

for every decomposition A^i + = -^i then its free energy density 
increases with N . We prove such a condition for a wide class of spin 
models which includes the Curie- Weiss model, its p-spin generaliza- 
tions (for both even and odd p), its random field version and also the 
finite pattern Hopfield model. For all these cases the existence of the 
thermodynamic limit by subadditivity and boundedness follows. 
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1 Introduction 



The rigorous theory of the thermodynamic hmit which aheady in the six- 
ties was a well established part of equilibrium statistical mechanics |R,uj 
received recently a new impulse thanks to the treatment of the Sherrington- 
Kirkpatrick |MPVj model of the mean field spin glass done by Guerra and 
Toninelli |GuToj . Moreover in a sequel work |Guj it became clear that a good 
control of the limit, especially when it is obtained by monotonicity through 
subadditivity arguments, may lead to sharp bounds for the model and carries 
important informations well beyond the existence of the limit itself. In this 
paper we build a theory of thermodynamic limit which apply to a family 
of cases including both random and non-random mean field models like the 
Curie Weiss model |Ba^, its p-spin generalizations, its random field version 
MPj . and also the finite pattern Hopfield model |Hoj . Due to the explicit 
size dependence of the local interactions we stress that mean field models do 
not fall into the class for which standard techinques [Ruj can be applyed to 
prove the existence of the thermodynamic limit. Moreover even in the Curie 
Weiss model in which the exact solution is available it is interesting to obtain 
the existence of the thermodynamic quantities without exploiting the exact 
solution (see [ENHnrTT] ). 

With respect to the theory relative to the random case jCDGGj we use 
here a different interpolation technique which works pointwise with respect 
to the disorder. The novelty of our approach relies on the fact that while 
in the previous case the condition for the existence of the limit is given 
in terms of a suitably deformed quenched measure, in the class of models 
we treat here we are able to give a condition with a direct thermodynamic 
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meaning: the Bolztmann-Gibbs state for a large system provides a good 
approximation for the subsystems. The fact that our existence condition is 
fully independent from the interpolation parameter relies on the convexity of 
the interpolating functional, a property still under investigation for the spin 
glass models ^nUEH^. 

2 Definitions and Results. 

We consider a system of N sites: {1, 2, A^}, to each site we associate a 
spin variable cxj taking values in {±1}. A spin configuration is specified 
by the sequence a = {ai,a2, ■ ■ ■ jCTn} and we denote the sets of all spin 
configurations by S^v = {±1}^. We will study models defined by a mean 
field Hamiltonian, i.e. for a given bounded function g : [—1, 1] —>■ M, 



Definition 1 For each N and a given inverse temperature (3 we introduce 
the partition function 




(1) 



where 




(2) 



j=i 




(3) 



the free energy density (and the auxiliary function a^) 

and, for a generic observable 0{a), the Boltzmann-Gibbs state 



(4) 





(5) 
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Remark 1 Obviously from ^ and one has that if g and g' are two 
functions from [—1,1] to M such that 

\\g-g'\\ := sup{|5((x) - g\x) \ : -1 < x < 1} (6) 

is hounded, then one has 

Wn - a^l < I3\\g - g'W (7) 

for every N . 

We can now state our main result: 

THEOREM 1 Let HN{a) be a mean field Hamiltonian (see eq. (OJj, l^^j- 

If for every 

partition of the set {1, 2, A^} into {1, 2, A^i} and {Ni + 1, A^} with 
N = N1 + N2 and 

Hni = HwiicTi, o'Ni) , = H^^i'^Ni+i, ■■■,o-n) , (8) 

the condition 

ujn{Hn) > ujn{Hn, + Hn,) , (9) 
is verified, then the thermodynamic limit exists in the sense: 

lim ajv = infaTv = «• (10) 

3 Proof 

Definition 2 Let us define the interpolating Hamiltonian as a function of 
the parameter t G [0, 1]; 

H^it) = tH^ + (1 - t)[H^, + H^,] , (11) 
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and consider its relative partition function Z]\f{t), free energy density f^it) 
and Boltzmann state uJN,t- 



The interpolation method that we are going to use is based on the sign control 
for both the first and second derivative of ajqit). More precisely the following 
holds: 

Lemma 1 Let H]\f be the mean field Hamiltonian and Hj\f{t) its relative 

interpolation. 

If 

jfN{t) < (12) 



for all t G [0, 1], then 



"Af < ^"A^i + ^"^^2 > (13) 



for each decomposition N = Ni + N2. 

Proof: trivially follows from the fundamental theorem of calculus and from 
the observation that definition (|TT|) implies: 

ZNil) = Z^ , (14) 

aAr(l) = On , (15) 
Z^iO) = Zn.Zn, , (16) 

and 

"w(0) = -^OiN^ + ^«^2 • (17) 

□ 



Lemma 2 Computing the t derivative of a]^{t), we get: 

= —-j^^N,t[HN — Hj^^—HN^]. (18) 



Lemma 3 The second derivative of a]^{t) is positive: 

= ^ajv(t) >0, (19) 
Proof: a direct computation gives 



= ^ {uJN,t [{Hn — Hn^ — Hn^Y~\ — ujIj^^ [Hn — Hn^ — Hn^]^ . (20) 

From Jensen's inequality applied to the convex function x i— > x^, it follows 
that a%{t) > 0. □ 

We are now able to prove the statement of Theorem 1. 

Proof of THEOREM 1. From Lemma © we notice that the hypothesis 

® 

is equivalent to the condition a^(l) < 0. On the other hand from Lemma 
it follows that a^(t) is an increasing function of t. This means that 
the determination of the sign of a'^{t) can be in general established by the 
evaluation of the sign in the extremes of the interval [0,1]. In particular we 
have: 

«^(1) < =^ a'^{t) < 0, Vt G [0, 1] (21) 



Using now Lemma (HJ, the subadditivity property fT^ holds for aN and 
then, by standard arguments |Ruj . 

hm aN = infaAT (22) 

Af^oo N 

The existence of thermodynamic hmit finally follows from boundedness of 
the function g in Eq. Indeed, calling K the maximum of g{x) on the 

interval [—1,1], we have 

aN = -\nY e^^'simM) > 1 ^PNk ^ _ /23) 

4 Applications 

In this Section we identify a class of mean field models for which the hy- 
potheses of our theorem are verified. Specifically these will be all models 
such that the function g of formula (0) is convex or polynomial. 

Corollary 1 Let the Hamiltonian he of the form 

Hn{a) = -Ng{mN) (24) 

with g : [—1, 1] ^ M a hounded convex function. Then the thermodynamic 
limit of the free energy exists. 

Proof: For a given cr G Sat and for every decomposition = A^^i + A'"2 we 
define the quantities 

m^,((T) = ^Xl^^ mN^a) = — ^ di , (25) 
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so that the total magnetization is a convex hnear combination of the two: 

^1 ^ ^2 .... 

Using this definitions the hypothesis Q of Theorem ((T]) is verified: 

f Ni N2 \ 

ujn {Hn - Hn^ - i^ivj = -Nujfq ( g{mfq) - —g{mfq^) - —giniN^) ) ^ 

(27) 

where the last inequality follows from convexity of g. □ 



Remark 2 The previous Corollary can be obviously generalized to the case 
where the function g is a convex bounded function of many variables, each 
of them fulfilling the property 

Corollary 2 Let the Hamiltonian be of the form 

H^ia) = -Ng{mM) , (28) 

with g : [—1, 1] ^ R a polynomial function of degree n G N. Then the 
thermodynamic limit exists. 

Proof: First of all we consider the case g{x) = (the generalization will 
follow in a simple way) with associated Hamiltonian 

N 



HNicr) = -Nm% 



fc-i 



«l,«2,---,«fc = l 



By splitting the summation into two pieces, the first containing the summa- 
tion with indexes all different among themselves, the second containing the 
remaining terms, we have 

1 



8 



(30) 



where the second summation J2* includes all terms with at least two equal 
indexes. A simple computation shows that 



E ^n---^^. = ^(1) (31) 

ii,...,ik 

Defining now the model with Hamiltonian 

1 ^ 
^^^^^ = -(iV-l)(iV-2)...(iV-A: + l). .E ^n^^.---^.' (32) 

it follows that 

H.W) - (iv-i)(iv-^);;-(iv-^.i) ^„(,) , 

= HN{a) + 0{l) (33) 

Using Remark 1 one has that the two models Hn and Hn have the same 
thermodynamic limit (if any). On the other hand for the model if^r we have 

= -(AT - l)(iV - 2). ..{N-k + 1) ^^^E^ J^(^^^^^^ ■ ■ -^^J 
= -NujNi(^icr2 - ■ -CTk) (34) 

where the last equality follows from permutation invariance ( ^^(ai-^o-i^ ■ ■ ■ ai^^) 
does not depend on the choice of the indexes). Analogously one can repeat 
the same computation for and in the state ujn- By permutation 
symmetry, this yields that Hypothesis Q of Theorem is verified as an 
equality 

ujN[HN-HN,-HN.)j=Q, (35) 

implying the existence of thermodynamic limit for the model iJjv, and so for 
the model if at. 
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Since we have proved the Corollary for g{x) = , the case of a generic 
polynomial function of degree n 

n 

(7(x) = J]afcx' (36) 

fc=0 

is treated by the same argument, applied to each monomial of the sum. □ 

Corollary 3 Using Remark 1 and the Stone- Weierstrass theorem, the ther- 
modynamic limit of exists if, instead of g being polynomial, g is merely 
continuous up to the boundary of [—1, 1]. 

4.1 Examples 

1. The Curie- Weiss models. 

For every integer p, with p < N, consider the model defined by 

1 ^ 

Hn{(t) := -j;^ ai^ai^---ai^ (37) 

11,22, .••,*p = l 

which represents the generalized Curie Weiss model with p-spin inter- 
action. The standard Curie- Weiss model corresponds to the case p = 2. 
From eq. (j2I) the previous Hamiltonian can be written as 

Hn = -NmFj^ (38) 

and the existence of the thermodynamic limit is then implied by Corol- 
lary Moreover, the same result holds for any linear combination 
(see Eq.(j2ni)) of generalized Curie Weiss models with p-spin interaction, 
both ferromagnetic and antiferromagnetic. 
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2. The random field Curie Weiss model. 

Here we consider the model defined by (see |MPj for a review) 

^ TV N 

HnQi, 0") := ^ o-iO-j + ^ hiai (39) 

i,j=l i=l 

where {/ii}i=i,...,Ar is a family of i.i.d. Bernoulli random variables, with 
probability distribution 

f 1/2, if hi = 1, 
Pih.) = { (40) 
[ 1/2, if hi = -1. 

For a given realization of the random field h, we define the quantities 



^U^^h) = ^J^^^-Y^''^ (41) 

i=l 



1 J^l-hi, 



m..(a,/i) = -^^-^ai (42) 

i=l 



The Hamiltonian can be written in terms of these variables as 



Hn = -Ng{m]^,m]^) (43) 

where 

g{m]^, m^) = (m^ + m^)^ - {m% - m^) (44) 

Since this function is obviously convex with respect to both m~l^ and 
m^, bounded by 2, and 



m^(a, h) = —m^^{a, h) + ^m^^(a, /i) (45) 



the Corollary ^ can be applied and we find (pointwise in the /I's) 

OiN{h) < -j^cuNiih) + —ttNiih) . (46) 
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Averaging now over the /I's, the subadditivity property for quenched 
aN is proved and this yields 



3. The Hopfield model. 

The Hamiltonian of the Hopfield model (see |Boj for a review) is given 
by : 

M N 

where M is the (fixed) number of pattern and the {Ci}i=i'.'.7N ^ 
family of i.i.d. Bernoulli variables with probability distribution 

f 1/2, if f = 1, 
= <^ (48) 
[ 1/2, ife = -l. 

Defining the quantities 

1 ^ 

<i^^0 = j;^J2^'''^ ,V/i=l,...,M, (49) 

1=1 

the Hamiltonian ()47|) can be written as 

M 

HM{a,0 = -Nj2«i^'0?- (50) 

The model can be included in the general treatment of the previous 
section by considering a function g oi M variables, 

M 

g{m],{a, 0, • • • , m^J{a, 0) = $^(m^^(a, 0)' (51) 

such that 

= -Ng{m],,...,m^). (52) 
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Since 

m%{a,0 = -^m%^{a,0 + ^rn%MO V/i = l,...,M (53) 

and the function g is convex with respect to every and bounded by 
M, using Corollary ^ we have (pointwise in the ^'s) 

«iv(0 < -^c^nAO + ^o^nAO ■ (54) 
Averaging over the ^'s yields the ()22|) . 

Remark 3 We want to notice that the method shown doesn't apply to 
the Hopfield model with a thermodynamically growing number of pat- 
terns defined for every positive constant 7 by 

-yN N 

Hn{^,^)--=-Y.nJ2^'^j^^^^^ (55) 
because this Hamiltonian is not of the form |2p. 
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